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Abstract 

We discuss the Penner type matrix model recently proposed by Dijkgraaf and Vafa for a possible 
explanation of the relation between four-dimensional gauge theory and Liouville theory by 
making use of the connection of the matrix model to two-dimensional CFT. We first consider 
the relation of gauge couplings defined in UV and IR regimes of Nf = 4, Af = 2 supersymmetric 
gauge theory being related as guv = "^(^hOV^^ir) 4 - We then use this relation to discuss the 
action of modular transformation on the matrix model and determine its spectral curve. 

We also discuss the decoupling of massive flavors from the Nf = 4 matrix model and derive 
matrix models describing asymptotically free M = 2 gauge theories. We find that the Penner 
type matrix theory reproduces correctly the standard results of M = 2 supersymmetric gauge 
theories. 



1 Introduction 



Recently an interesting observation has been made on the relationship between the M = 2 
supersymmetric gauge theory and the 2d conformal field theory. In [T], it has been shown 
that Nekrasov's partition function [2] of M = 2 SU(2) superconformal gauge theory can be 
identified with the correlation function of 2d Liouville theory. This relation was further studied 
and checked from various directions [3j H]. For general discussions on S-duality in M = 2 
superconformal gauge theories, see [51 El [7]. 

In ref. [8], the following matrix model with a Penner like action 

3 

W(M) = J2^ lo S( M " ?! = 0, q 2 = 1, q 3 = q (1.1) 

i=l 

has been proposed to describe M = 2 SU(2) superconformal gauge theory with four hypermul- 
tiplets. It is suggested that this matrix theory explains the correspondence of pQ, by making 
use of the CFT description of the matrix model 0, EO]- (This matrix model corresponds to 
c = 1 (or b = i) case in the Liouville theory.) The four mass parameters of the gauge theory 
are identified with /Xj in the matrix model action and /io which corresponds to the charge at 
infinity. These parameters satisfy the following relation [8] 

3 

^2(* + IM> = -2g a N, (1.2) 

8=1 

where N is the size of the matrix. 

We study this matrix model in details in this article and show that they in fact correctly 
reproduce the physics of SU (2) Seiberg-Witten theory [HJ [12] . 

In section [2] we first consider the M-theory curve for the Nf = 4, M = 2 supersymmetric 
SU(2) gauge theory based on the brane construction [13]. We study its period integrals and 
by comparing with those of the standard Seiberg-Witten curve we obtain the relation between 
gauge coupling constants at UV regime and IR regime as guv = ^(^hOV^^ir) 4 - The pa- 
rameter q which appears in the matrix model action (11.11) is identified as the UV coupling 

q = guv- 

We discuss in section [3] the modular property of the matrix model (11.11) by making use of 
the IR-UV relation of gauge coupling constants. We study the modular transformation of the 
spectral curve and determine its precise mass- dependence by imposing modular invariance. 

In sections HI we consider the decoupling of massive flavors and find matrix model actions 
corresponding to asymptotically free gauge theories with Nf = 2, 3. We check that the discrim- 
inants of the spectral curve and free energy of these matrix models match well with those of 
Seiberg-Witten theories. 
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2 Af = 2 gauge theory and UV and IR gauge coupling 
constants 



2.1 M-theory curve of Nf = 4, SU(2) gauge theory 

M = 2 SU (2) gauge theory with four hypermultiplets is known to be scale invariant and has 
the exactly marginal coupling constant 

Tuv = \~ -o— • (2.1) 

The flavor symmetry of the theory is 50(8) whose maximal subgroup is 5f/(2) 4 . While in- 
troducing the masses of the hypermultiplets breaks the conformal invariance, the theory is 
modular invariant involving the triality of 50(8) which rotates the mass parameters [12]. We 
introduce the mass parameters of four hypermultiplets m± and rh±. By combining these, we 
define 

m = ^(m + -m_), m 2 = ^(m + + m_), m 1 = ^(rh + -m_), m 3 = ^(m + + m_), (2.2) 

each of which is a mass parameter associated with each SU (2) flavor symmetry. 

This gauge theory can be obtained by considering the intersecting brane system in type IIA 
string theory [13] • The SU(2) gauge part is induced by two D4-branes suspended between two 
NS5-branes. The D4-branes occupy the x 0,1,2,3 and x 6 directions and the NS5-branes occupy 
2.0,1,2,3 anc j x 4,5 are denoted by the complex coordinate v. A massive hypermultiplet 

can be introduced by a semi-infinite D4-brane attached to an NS5-brane. We here choose the 
configuration such that two semi-infinite D4-branes are attached to the left NS5-brane and 
extended to x 6 = — oo and two more D4-branes are attached to the right NS5-brane extending 
to x 6 = +oo. By the M-theory uplift, the Seiberg-Witten curve [11] of this theory becomes a 
hypersurface on (t,v) G C 2 where t = e -( x6 + lxW )/ R : 

(v - m+)(v - m_)t 2 + d(v 2 + Mv - U)t + c 2 (v - rh+)(v - m_) = 0, (2.3) 

where M and U are constants which depend on the masses and the Coulomb moduli u. The 
first and the third terms are determined as follows: in the large t, the first term is dominant 
and we obtain v ~ m+, m_ which should be the masses of the hypermultiplets induced by the 
left semi-infinite D4-branes (large t corresponds to x 6 —> — oo). On the other hand, in the small 
t, the third term is dominant and v ~ m + ,m_ which are the mass parameters induced by the 
right semi- infinite D4-branes (small t corresponds to x 6 — > +oo). The dimensional analysis and 
the regularity constraint in the massless limit show that M is linear in the mass parameters 
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and U is linear in the Coulomb moduli parameter u and also contains additional terms in mass 
squared. Also, c\ and C2 are the constants which parametrize the gauge coupling constant. 
Then, the curve can be written as 

v 2 {t - l)(t - q) = v((m+ + mJ)t 2 + (1 + q)Mt + g(m+ + in J)) 

—m + m_t 2 — (l + q)Ut — qm + m~, (2-4) 

where we have chosen that c\ = — (1 + q) and C2 = q. 

By eliminating the terms linear in v and changing the coordinate as v = xt, the curve can 
be written as the following form [5]: 

2 f m 2 t 2 + (l + q)ft + m 3 q \ 2 (m 2 - m 2 )t 2 - (1 + q)Ut + (mj - m 2 )g 

where a; 2 <it 2 is considered as a quadratic differential on a sphere (t is a coordinate on the sphere) 
and has double poles at t — 0, 1, q, oo. In this coordinate, the Seiberg-Witten one-form is 

Asw = u^c (2 - 6) 

where we have divided by the factor 2\/27ri in order to be consistent with the convention in 
later sections. Here, (t, x) are local coordinates on the cotangent bundle of the sphere. The 
Seiberg-Witten curve is the double cover of this sphere with four punctures. 

The moduli space of the sphere with four punctures is parametrized by q above. As discussed 
in [T3| 15], this moduli space can be identified with the parameter space of the exactly marginal 
operator of 4d gauge theory, that is the gauge coupling constant (12. lft . Therefore, we identify 
q with guv = e nlTUV . Note that this is the UV gauge coupling constant. The effective gauge 
coupling at IR is derived using the full information of the Seiberg-Witten theory. We will 
discuss this point in the next subsection. 

The masses of the flavors can be read from the residues of Asw at the punctures. As seen 
in [5] , each puncture is associated with each SU (2) flavor symmetry. Indeed, we can easily see 
that the residues at t = 0, oo are ±mi and ±mo where the signs represent their value on the 
upper and the lower sheets respectively. (More precisely, the residues are ± 2 ' we w ^ 
ignore the factor 2\^2iri in what follows when discussing the residues.) Also, the requirement 
that the residues at t = 1, guv should be ±7712, ±1713 determines the form of M as follows: 

M = _ 2guv_ (2 7) 

1 + guv 

The Coulomb branch vev a can be obtained by the period integral 

Asw- (2.8) 

A 



This is the function of the Coulomb moduli u. Also, the dual B cycle integral computes the 
derivative of the prepotential 

dJ 7 f 

a D = -q^ = J A sw - (2.9) 

The moduli derivative of Asw is a holomorphic one-form 

dA S w (1 + guv) dU 

oj = — — = . — — dt, (2.10) 

du 4m^/2P^t)du v ; 



whose A and B periods give |^ and P±(t) is a degree- four polynomial whose precise form 
does not concern us. As we will see in subsection 12. 2\ ^ becomes 1 in the weak coupling limit: 
guv -> 0. 

We now consider a way to extract the parameter U (essentially (Tr<j) 2 )) from the M-theory 
curve above. In the pure M = 2 SU(N) gauge theory, it is known that such an operation exists 



(Ti<P 2 ) = ^-^v\ sw = ^-^x 2 tdt, (2.11) 

where v is the coordinate in (12. 3p . Also in the SU(2) gauge theory with four massive flavors, 
we can use this kind of integral to extract the U parameter. Actually, since the term including 
U in x 2 is 

0- + quv)U 



i(t -l)(*-<Zuv) 

the integral 



(2.12) 



— & x 2 t n dt, for n > 2 (2.13) 



2m 



C ry 



is linear in U and contains terms involving mass parameters. We denote a contour around 
infinity in the counterclockwise direction as C^. In the simplest case (n = 2), we easily obtain 

— 4 x 2 t 2 dt = (l + q VY )U-{l + q VY )m 2 + {q V y-l)m 2 2 + 2q V ym 2 m 3 . (2.14) 

This relation becomes an important point when we compare the gauge theory with the matrix 
model in section [3j 



2.2 UV and IR gauge coupling constants 

The M-theory curve (12. 5p which we have seen in the previous subsection is different from the 
standard Seiberg-Witten curve introduced in [T2]. In this subsection, we will see that the 
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comparison of the M-theory curve and the standard Seiberg-Witten curve leads to a relation 
between the UV and IR gauge coupling constants [HJ [THl [TJ H]. We first consider the massless 
case for simplicity. 

On the one hand, we have the standard form of the Seiberg-Witten curve [12] 



y 2 = 4x 3 - g 2 u 2 x - g 3 u 3 , (2.15) 

where u parametrizes the Coulomb moduli and is related to (Tr0 2 ). The IR gauge coupling 
constant 

_ 6 m 8iri 

tir = h — (2.16) 

is calculated from the period integrals of the holomorphic one-form u which is defined as 

uj = = (2-17) 

47T y Air ^4 X 3 - g 2 u 2 x - g 3 u 3 

where g 2 and g% are the functions of giR = e mnn : 

92(e>i, qm) = ^ (MQm) 8 + Mlmf + M<lm) 8 ) , 

03 (wi, qm) ={£) ^ (Mqm) 4 ~ Mool)*) ^(gm) 8 + M^Mool)*) , ( 2 -18) 



and 2ui = ir in this case. A-cycle integral of u gives 



da y/2 f dx 1 

du Air J A ^x 3 - g 2 u 2 x - g 3 u 3 2y/2u 



(2.19) 



which leads to a = hv2u. 

As pointed out in [T7|, however, the gauge coupling constant receives the correction due to 
instanton effects even in Nf = A theory and, therefore, giR ^ guv- As we saw in the previous 
subsection, we obtain the following curve in the massless limit of the M-theory curve: 

t(t-l){t-q vv y v ; 

where guv is the UV gauge coupling constant. Note that the parameter U in the above curve 
would be different from u in the curve ( 12. 15ft . Dimensional analysis shows that U is proportional 
to u: 

U = Au, (2.21) 
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where A depends only on the UV gauge coupling guv- The holomorphic one- form is given by 



w = i >a+^ * (2 . 22) 

AV2m V ix v/t(t - l)(t - g uv ) 

By changing the coordinate as t — > t + 1+ g UV and by rescaling as t — > 42, £(£ — — g uv ) is 
transformed to 

16 ^ 3 - ^(1 - g uv + g uv )^ - ^(2 - 3guv - 3g uv + 2g uv )) = l%o- (2-23) 
Then, the holomorphic one-form reads as 



u 1 /(H^, ( , 24) 
4V27T V it 2/0 

Now, we want to compare this with the one- form (12.171) . In order to reproduce the A period 



— 5=, we should obtain 

2v Zu 



dz 2n 



T A y y/(l + g uv ) A 
Therefore, the following relations must be satisfied: 

Y^C 1 -?uv + 9uv) =02(^1, flu), 



This shows that 

Ami = 

#3 (?m) 

and 



(2.25) 



-^(2 - 3g uv - 3g uv + 2g uv ) = 93(01, ?m), (2.26) 
and the period is given by 

2wi = ; 2n (2.27) 
v/(l + guv)^ 

By using formula $\ = $| — -d^ the right hand sides of (I2.26|) can be calculated as 

S2= ^ (1 + 3Dv) ^i( 1 -| + (| ^ 

*^(l + ? uv)^(2-3j-3(J) 2 + 2(jV|- ("8) 



9uv = ! 2 i 9lR |! = 16gm - 128 ft 2 R + 704g I 3 R - 3072g I 4 R + . . . . (2.29) 



A ~ Mmylumr (2 ' 30) 
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which implies that U = u + 0(q\jy) = (Tr 2 ) + C(guv) m the weak coupling limit. The relation 
(12.291) between the UV and IR gauge coupling constants has already been suggested in [18] from 
the topological string analysis and derived in [I] along the similar line as above. 

Let us next consider the case with the massive flavors. In this case, the parameter U 
appeared in the M-theory curve ( 12.5ft receives additional contributions due to mass terms. 
However, the form of A (12.301) should not change. Therefore, we obtain 

dU 1 



It follows from this that the holomorphic one-form ( 12.101) becomes 



(2.31) 



dt 

u = ; 2.32 

The A-cycle integral of this one-form can be translated into the following standard form of the 
elliptic integral of the first kind, by the coordinate transformation: 

du 

(2.33) 



y/(l-U 2 )(l-k 2 U*y 

up to a constant factor, where k is a function determined from P±(t). We also obtain the similar 
form for the B cycle integral. The IR gauge coupling constant is given by the period of this 
one form: Tra = l 2 -^ and this in turn leads to 

This determines the IR gauge coupling constant in terms of the UV one. 



3 Matrix model and modular invariance 



In this section, we analyze the matrix model with the action (11.11) which was proposed to 
describe J\f = 2 SU(2) gauge theory with four massive hypermultiplets and study its modular 
properties. We first review the general technique of the one matrix model and then consider 
the modular invariance of the spectral curve. 

We define the free energy of the matrix model as follows: 



exp 



9 2 s 



{ — W{M) ) 
9s J 




exp 



XV(Ai) + <7.X>e( A '- A - 



KJ 



(3.1; 
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where M is the hermitian matrix whose size is N and W(M) is (11. ip . In the second line, we 
have switched to integrals over the eigenvalues A/. The last term comes from the Vandermonde 
determinant. The critical points are determined by the equation of motion 

By ignoring the second term, we obtain two critical points e p (p = 1,2). The positions g» will 
be chosen as q± = 0, q-i = 1 and g3 = q. Let each N p (p = 1, 2) be the number of the eigenvalues 
which are classically at e p . Let us define the resolvent of the matrix model as 

R m {z) = g s tr ^ . (3.3) 

In what follows, we take the large N limit while \x v and the filling fractions = g s N p are 
fixed. The loop equation in the large N limit is written as 

(R m (z)) 2 = -(R m (z))W'(z) + M (3.4) 

where 

W'(z) - W'(M) \ 



has simple poles at z — qi. Note that this function is not a polynomial in contrast to the 
ordinary case of polynomial action, q are functions of /ij and the filling fractions but they 
satisfy ^E c, = which follows from the equations of motion (13. 2p . 
Let us define the meromorphic one-form A m = ^ z l dz . such that 

x[zf = (2(R m (z)) + W'{z)) 2 = W'{zf + f(z). (3.6) 

This one-form has simple poles at z — q^oo with the residues Hi, fiQ. (As in the field theory 
analysis, the residues are precisely - by the definition of X m . But, for convenience, we will 
ignore the factor in the denominator when discussing the residues.) Note that the residue of 
the pole at z = oo can be evaluated by observing (R m ) ~ and W'(z) ~ — at large z 
and by using the relation (jl.2p . 

The spectral curve of the matrix model (I3.6P looks like the M-theory curve. The former 
includes the parameters fii which will be identified with the masses m ; of the latter. These 
are the residues of the meromorphic one-form X m and Asw & t the simple poles. Hence, the 
identification q = guv is needed. To compare more precisely, we should add a dimensionful 
parameter by g s — > g s e since the parameters in the matrix model are dimensionless. In the 
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analysis of the gauge theory, the parameter U is obtained from § x 2 t n dt for n > 2 (I2.14p . In 
the matrix model language, the similar integral gives 

<f x 2 z 2 dz — ~~. & (2R m + W'fz 2 dz = 4/i ^(trM) + 2/i (/i 2 + ^ 3 ), (3.7) 

where is the contour around infinity. Therefore, the precise correspondence between the 
matrix model and the gauge theory is given by the following identification: 

Afi g s (trM) + 2^{^ 2 + qnz) 

= (1 + qvv)U - (1 + qvv)ml + (<?uv - 1)m\ + 2q vy m 2 m 3 , (3.8) 

and also ^ — rrii, q — guv- 

We can rewrite (trM) in terms of the coefficients q in (I3.5p . First of all, the residue at 
z = oo imposes the following constraint on q: 



c 2 + qc 3 = l4 



X> ■ (3-9) 



j=i 



Therefore, by recalling £\ c i = 0, onr y one °f c / s is independent. By using (13 .9p and q = 0, 
the spectral curve can be written as 

j = ((± + J^ + J*Y + ^o-fi^)^ + ^, (3.10) 
\z z—1 z — qj z(z — l)(z — q) 

Also, the filling fractions can be obtained by the A-cycle integrals of the one-form: 

A m , p=l,2, (3.11) 



where A p are the cycles around the branch cuts (corresponding to two critical points). From 
this, we can in principle determine v p as the function of /i, and c\. In other words, c\ can be 
treated as an independent parameter. 

It follows from the dimensional analysis that c\ has mass dimension two. Therefore, it 
is natural to consider that c\ is a linear function in U of the M-theory curve. Indeed, by 
substituting the explicit expression (13.10p into the left hand side of (13.711 . we obtain 

— * x 2 z 2 dz = -qa - (1 + q)nl + (1 + q)& + (q - 1)^1 
21X1 Jc^ 

+ (1 - q)(4 + 2qn 1 [i 2 + 2/x 3 /ii. (3.12) 

This and (13.81) lead to 

?uvci = (1 + qv\)m\ + (1 - <?uv)™3 + 2q VY m 1 m 2 

— 2q VY m 2 m 3 + 2m 3 m 1 - (1 + quv)U. (3.13) 
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We can show that under this relation, /ij = and q = guv, the meromorphic one- form A m is 
equal to Asw- Therefore, the A and 77-periods of A m and Xsw coincide, e.g., \/2vi(u) = a(u). 
Also, w-derivative of A m gives a holomorphic one-form: 

= ^ = ^ (3.14) 

where Pm^z) is the same degree 4 polynomial as P 4 in ( I2.10p . The effective gauge coupling 
constant tjr can be obtained from the period of oj m . 

3.1 Modular invar iance 

It is known that the standard Seiberg-Witten curve of 577(2) gauge theory with Nf = 4 is 
invariant under modular transformation in [12]. We will see in this subsection that the spectral 
curve (13.101) which was identified with the M-theory curve (I2.5P can be made modular invariant. 
To begin with, we consider the massless limit of the curve (I3.10j) 

u 

o_ (1 + gUvg _ g m ,, 

z(z-l)(z-q)- z(z-l)(z-qY ^ 0> 

where we have used the relation f)2.30p . This is invariant under the following transformations 

I : (z, x) ->• (1 - z,x), g uv ->■ 1 - guv, u -u, S (3.16) 

II : (z,x) ^ (-,-z 2 x), guv -> — , u->u, STS (3.17) 
z guv 

Since guv = ^2/^3 as seen m subsection [221 the former transformation is the exchange of $2 an d 
t?4, which is the S-transformation. Also, the latter one can be seen as the STS-transformation: 
■&\ <-> We note that the sign of u changes under S-transformation. This behavior is the 
same as in the case of standard SW curve. In fact in ( I2.15P and g% are even and odd and u 
changes sign under S-transformation. 

Next, let us consider the massive case. As analyzed in [12], the modular transformation in 
this case involves the triality of 50(8) which permutes the SU(2) flavor symmetries. In our 
notation, this permutes the mass parameters m, (i = 0, . . . , 3) associated with these 5f7(2)'s. 
Under the S- and STS-transformations, positions and residues of the poles of A m are transformed 
as 

7 : (0,1, g uv , 00) (1,0,1 -g uv , 00), m 1 ^ m 2 , (3.18) 
77 : (0, 1, g uv , 00) ->• (00, 1, ,0), m ^ m 1 . (3.19) 

guv 
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Under these transformations, the spectral curve should be invariant. Let us substitute ( 13.1 3p 
into ( 13.101) and we obtain 

2/ v f m i m 2 m 3 \ 2 (mjj - (Eti m) 2 )z + (1 + guv)m 2 
x (z;mi;q vv ) = 1 H H (3.20) 



2 z — 1 2; — guv , 

+ (1 - guv)"i 3 + 2quym 1 m 2 - 2qu V m 2 m 3 + 2mim 3 - (1 + quv)U 

z{z -l)(z- guv) 

We then impose the conditions 

£ 2 0; m , m 1 ,m 2 , m 3 ; guv) = ^ 2 (1 - 2; m , m 2 , m u m 3 ; 1 - g uv ), (3.21) 

x 2 (z;m ,m 1 ,m 2 ,m 3 ; g uv ) = -jX 2 (-;m 1 ,m ,m2,m 3 : ). (3.22) 

z guv 

Requirement of modular invariance determines completely the mass dependence of the param- 
eter U. The solution to the above conditions is given by 

(1 + Quv)U = ^ - guvK + m 3 ) 2 + (3 23) 

3.2 Relation for (trM m ) 

In (I3.8p . we have written down the parameter U in terms of (trM) by using the integral § x 2 z 2 dz. 
As seen in subsection 12.11 we can in principle use the other integrals as 

— / x 2 z n dz, for n > 2, (3.24) 

in order to extract the parameter U . This in turn gives the relations for (trM m )'s. For example, 
we consider n = 3 case. From the M-theory curve, we can compute 
1 



x t dt = 4m (l + q vw )g s {tiM) + guv^o ~ <?uv™i + (guv - l)m 2 

2m J r 

+guv(l - guv)™ 3 + 2(1 + q\jv)m (m 2 + guv«i 3 ), (3.25) 

where we have used (13. 8p . On the other hand, the same integral for the matrix model is 
calculated as 

1 



2m 



x 2 z 3 dz = 4^ g s (tiM 2 } - 4g s (trM)(g s (trM) + m 2 + q vy m 3 ) 

Coo 

-{m 2 + guv™ 3 ) 2 + 2m (m 2 + gu V m 3 ). (3.26) 

From these equations, we obtain 

Afi g s (tiM 2 ) =Ag s (tiM){g s (tiM) + (1 + g UV )™o + m 2 + q vy m 3 ) 

+quv(ml -m\ + m\ + m 2 ) + 2q\jy(m 2 m 3 + m m 2 + m m 3 ). (3.27) 

Similarly, (trM m ) for m > 3 can be written in terms of the lower order ones. 
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3.3 Free energy 



So far, we have learned that M-theory curve and spectral curve of matrix model can be iden- 
tified. However, it is not so obvious that the matrix model free energy coincides with the 
prepotential of gauge theory, in other words, the B-period of A m (differential of the matrix 
theory) is written as Let us see this below. We consider the free energy and its derivative 
with respect to the filling fractions. First of all, we rewrite the free energy as 



e F m /gl _ exp 



-^U d\p(X)W(X) + J d\dX'p(X)p(X') log |A - A' 
where p(X) is the eigenvalue distribution function normalized as 



(3.28) 



dXp(X) = g s N, 

or p(A) = g s ^2 j 5(X — Xj). The resolvent can be written as 

P(A) 



(3.29) 



Rm(z) 



dX 



z-X' 



(3.30) 



Now, let us consider the derivative of the free energy with respect to the filling fractions 
Up. Here, we will follow the discussion in JT9J [20]. The variations with respect to v v can be 
considered as the shift p(A) — > p(A) + (Su p )5(X — e+) where e+ are any points on the branch 
cuts (except for the end points of the cuts). Therefore, the free energy is shifted as 



SF m = (5u p ) 
= Op) 



H/(e+) + 2 J dXp(X)\og\X-e; 

W(e+) + ^-J (R m (z)) log(z - e+)dz 
im Ja 1 +a 2 



(3.31) 



The second term is evaluated by deforming the contour to the cycle around the log cut as 



2ni 



— i (R m (z))\og(z-e+) 



Ai+A 2 



(R m {z))dz + -— (b £- dz 

2m z 



{R m (z))dz + 9s N "(log A + m) + O Ij- . (3.32) 

Note that (R m ) has a pole only at z = 00. We regularized the integral by introducing the cut 
off Aq. We will take Aq to infinity after the calculation. By changing (Rm) into A m we obtain 



f A ° ( 1 

A m -/i logA + 27r^N + e> f — 



(3.33) 
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In order to compare with the prepotential of the gauge theory, we let g s — > g s e as above, 
which gives /ij, /iq, v p and F m dimension. Since F m has mass dimension 2, it is written as 



2F 



i=l 



p=l,2 



<9z/„ 



(3.34) 



The first term is evaluated as 



dF 



: 5s (Tr log(M - qi )) = ^-.l (Rm(z)) \og(z - q t 

2m J Al+A2 



)dz 



Til 

7! 



+ 



Ao 



Ao 



A,, 



-w( qi 



//() 



logA + g 3 Nm + O 



1 

Ao 



(3.35) 



where we have deformed the contour to the cycle around the log cut as in f 1 3 . 3 2 j) . Also, by using 
(13.331) . the second term in (I3.34[) is 



p=l,2 



' dVr, 



-V2 



Til 



-U 2 



Li, 



\ rn - g s NLi log A + 27ii(g s N) 2 + 0[ — 



(3.36) 



\ Y1 V* + ^° ] / Am + y A ™ - lo S A o + 2TTi(g s N) 2 + O 



where Bi and -B2 are the paths from A to A through the branch cuts respectively and we 
define the cycle B — B% — B\. Combining these, we finally obtain 



2F„ 



711 

7! 



i=l 



+V0 



B 2 



X m + ni (b A, 



A, 



B 



V* W &) + \& lo S A o " ™9sNfio + O (J^ . (3.37) 



The first line of the above equation is the same as the prepotential of the gauge theory 
up to a factor iri which can be absorbed by the redefinition of the free energy. Note that the 
divergence of integrals in the first line is canceled by the divergence of terms in the second line 
and the free energy is finite. There are also finite terms in the second line which depend only 
on the mass parameter \x% and liq. However, we should note that the terms in the first line 
have ambiguities due to the choice of integration paths. By deforming the paths, an additional 
contribution which is bilinear in the masses can appear. Therefore, we conclude that the free 
energy of the matrix model is the same as the prepotential up to these moduli independent 
terms. 
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4 Matrix model for asymptotically free theories 



4.1 Asymptotically free theory with Nf = 3 

In this and next subsections, we consider the decoupling of heavy flavors from the matrix model 
(1.1) and introduce the matrix theory with flavors Nf < 4 which would describe the M = 2 
gauge theories in the asymptotically free region. 

Before taking the decoupling limit we should recall the precise relationship between u and 
(Tr (ft 2 ) as pointed out in ref [12] 

u={Kcf)- l -{d\ + $t)j^ml (4.1) 

i=0 

By substituting (14.11) into (I3.23j) . we obtain 

(Tr 2 ) 1 1 

- (1 + quy)U = h 2q vv m 2 m 3 + -quv(m 2 3 - ml) + -q V \{m 2 2 ~ m l)- ( 4 - 2 ) 

In the following subsections we again use the symbol u to denote (Tr</> 2 ) for simplicity. We 
hope that no confusion should arise from this change of notation. 

In order to discuss the decoupling limit we write the mass parameters as follows 

m± = m 2 ± m , m± = m 3 ± m l5 (4-3) 

and consider the limit rh_ — > oo, guv — ^ and the other masses and m_guv = ^-3 fixed, where 
A 3 corresponds to the dynamical scale of the gauge theory. In this case, the relation on the 
mass parameters ( 11. 21) becomes 

m + m 2 + m+ = -2g s N, (4.4) 
and this remains finite. The matrix model action leads to 

W(M) =m + \ogM - ^ + m 2 \og(M -1). (4.5) 

and we obtain the spectral curve for Nf = 3 theory 

2 _ A| m+A 3 u-(m 2 + |m+)A 3 m\ m 2 A 3 , , 

X ~ Az 4 z 3 (z-l) z 2 (z-l) z(z-l) z(z-l) 2 z 2 (z-l)' [ ' } 

Note that the structure of the singularities is similar to that obtained in [B]. Let us consider 
M-derivative of the Seiberg-Witten one-form A r 



xdz . 



u _ dz -1 _ 1 

du 2y/2Tii2xz 2 {z - 1) ~ Am^2Q A (zy 
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where 

Q±(z) =mlz 4 + f-u - m 2 , + m 2 , + ^m + A 3 ^ z 3 

The discriminant of the polynomial above completely matches with that of the Seiberg-Witten 
curve for Nf = 3 in |12] : 

y 2 = x 2 (x — u) — -A 2 (x — u) 2 — ^{ m + + m2 - + ?7i+)A 3 (:r — u) 

+m + m_m + A 3 x — -(m+m 2 ^ + m 2 + m 2 + m 2 + )k\. (4-9) 

Note that the dynamical scale A 3 differs by a factor from Asw that appeared in [12]. This is 
due to the difference of the decoupling limit of ours and that of [12] where they fixed 

64m„g IR = A sw . (4.10) 

As discussed above, guv = 16gm for the weak coupling region. Thus, we obtain Asw = 
4m_g uv = 4A 3 . 

For the equal hypermultiplet mass case where m = 0, m 2 = m and rh + = m, the discrimi- 
nant becomes 

A 2 / 1 \ 3 

A = — - I u — m 2 — -mA 3 (-32m 2 + 2A 2 w - 48wmA 3 + 3mA| + 6m 2 A 3 + 128m 3 A 3 ) . 
32 V 2 / 

(4.11) 

This implies that three (SU (3) triplet) massless particles appear at the singularity u = m 2 + 
|mA 3 . This is the correct property of the curve (j4.9p . 

Free energy 

Now we consider the free energy of this matrix model. The formula (13.331) for the derivative of 
the free energy with respect to the filling fractions can be used without change. By restoring 
the dimension, we obtain the following equation 

2F - = y n -ik + m *ik + ^ + F ~ (412) 

\ P — / 

By using the result in the previous section, the computation of the first two terms are straight- 
forward: e.g., 

dF m in ( f° f A °\ 1 l T „ /ft x mo, » . Ar / 1 



dm + y/2 



J^+j_ °\ \ m -l W (0)-^\ogA + nig s N + O (J^J . (4.13) 
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Also, the last term is exactly the same as ( I3.36P except that Y2i=i m * becomes m + + m 2 . Only 
the new ingredient which need a little bit care is the third term in ( 14.12ft . In order to compute 
this, we note that in the spectral curve: x 2 = W'(z) 2 + f(z), f(z) is written as 

This /3 can be directly evaluated from the curve (14. 6p . which leads to 
dF m 9s / i\ f 3 1 / I 



<9A, 



f(Ei) = ^ = ^(« + (r + + ^)A3-^). (4.15) 



By collecting these, we finally obtain the expression for the free energy 



»0 /•! 



2F m = ym + J +m 2 / +m / )\ m + iri(j) X m (f> X m + " + . . . . (4.16) 



1 JB 2 / J Ax JB 4 



This nicely matches with the prepotential of the theory with Nf = 3 [2T1 1221 [23] up to the 
irrelevant factor iri. The ellipsis is the terms which depend only on the masses and the scale 
factor. It is reassuring to note that the beta function coefficient b = 1 in front of u is correctly 
recovered in the matrix model computation. 

4.2 Asymptotically free theory with Nf = 2 

In this subsection, we propose a matrix model action which corresponds to SU(2) gauge theory 
with Nf = 2. Method of its derivation is similar as the one in the previous section: we consider 
the limit where z is rescaled as j^z and m_ — > oo while keeping m„A3 = Af. The spectral 
curve is given by 

n Ni fh+Ao u m+Ao An 
Az 4 z 3 z 2 z 4 

Here, we have rescaled x as x — > ^-x. This curve is the same form as the first realization of 
SU(2) gauge theory with Nf = 2 in [6]. In this limit, the action of the matrix model (14. 5ft 
becomes 

W(M) = m + logM - ^ - (4.18) 

where we have ignored irrelevant constant terms. The relation on the mass parameters ( 14.41) 
is, in this limit, 

m+ + m+ = -2g s N. (4.19) 
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The holomorphic one-form is 



d xdz dz 
u = t; =— = — 7= ; (4.20) 



where 



A 2 A 2 
R A ( Z ) = -2-z 4 + m + A 2 z 5 + uz 2 + m + A 2 z + —2. (4.21) 

The discriminant of the above polynomial agrees completely with that of the Seiberg-Witten 
curve for Nf = 2 in |12j : 

y 2 = {x 2 - l -At){x -u) + m + m + A 2 2 x -\{m 2 + + rh 2 + )A A 2 . (4.22) 

For the equal mass case where m + = m + = m, the discriminant becomes 

A = A\ (u + 2mA Nf=2 + ^A 2 ,^ (u - 2mA 2 + (u - m 2 - ^Aj}j . (4.23) 

Free energy 

Let us compute the free energy of this model. The procedure is the same as that in section [3] 
and 14.11 Only non-trivial point is the calculation of A 2 derivative: 

£~f(i: -£-!<!>>■ (4 - 24) 

where g 2 is the coefficient of g(z) = — + % in the spectral curve: x 2 = W'(z) 2 + g(z). By 
comparing this form of the spectral curve with (14.171) . we obtain g 2 = u + |A 2 — rh 2 + . The 
second term is A/} = (trM). This can be determined, by using the same argument as that 
in subsection I3.2[ as 

2A 2 (trM) = — u — -A\ + m\. (4.25) 

Therefore, we obtain 

dF m If 1,, 9 \ u mh + mi - A\ . 



2F m = — ( m+ / +m+ / } A m + m <p X m <p \ m + — + ... . (4.27) 



Then, the similar computation as above leads to 

A m + 

'B 2 / ^Ai JB 

This is the same form as the prepotential of the theory with Nf = 2. It is nice to find the 
appearance of the beta function coefficient 6 = 2 for Nf = 2 theory in front of u from the 
matrix model computation. 
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5 Conclusion and discussion 



In this article we have studied the Penner type matrix model and gave strong evidence that 
it reproduces correctly the physics of M =2 supersymmetric SU (2) gauge theories. It will be 
interesting to see if the model can in fact be used to provide a simple derivation of the AGT 
relation. 

Generalization to SU(2) linear quiver gauge theories [5] may be straightforward. In this 
case, the corresponding action of the matrix model is 



where k corresponds to the number of the simple poles of the Seiberg-Witten one-form. It will 
be interesting to study the S-duality group of this matrix theory. 

A different type of matrix model has been proposed in ref [23] to derive Seiberg-Witten 
theory. It is important to see how this model is related to the matrix theory discussed in this 
paper. 
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